Purity as a witness for initial system-environment correlations 

in open-system dynamics 
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We study the dynamics of a two-level atom interacting with a Lorentzian structured reservoir 
considering initial system-environment correlations. It is shown that under strong system-reservoir 
coupling the dynamics of purity can determine whether there are initial correlations between system 
and environment. Moreover, we investigate the interaction of two two-level atoms with the same 
reservoir. In this case, we show that besides determining if there are initial system-environment 
correlations, the dynamics of the purity of the atomic system allows the identification of the distinct 
correlated initial states. In addition, the dynamics of quantum and classical correlations is analyzed. 
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I. INTRODUCTION 

The preparation of the initial state plays a key role 
in the implementation of any experiment. As recently 
noted, 1] different procedures of the preparation of ini- 
tial states can lead to nontrivial differences in the ex- 
perimental results due to interaction with the environ- 
ment. These differences arise from initial correlations 
between the system and its environment. Although the 
assumption of uncorrelated system-reservoir initial states 
is widely used, it is not always well justified, specially 
when the system strongly interacts with the environment 
Thus, the influence of initial correlations on the dy- 
namics of open quantum systems has been intensively 
studied [3l-[Tl|. In Ref. [l^], the authors verified that the 
linear transient response of a two-level system coupled to 
a reservoir differs significantly if the initial state is cor- 
related or not. Another important issue concerns about 
the possibility of associating a completely positive map to 
a quantum dynamics, which is possible only when there 
are no initial quantum correlations between system and 
reservoir [11]. Moreover, the effects of the initial correla- 
tions were recently observed experimentally by using an 
all-optical apparatus to probe the evolution of the two- 
qubit polarization entangled state [l^ . 

The usual way of measuring correlations in quantum 
systems is by mutual information ^is']. Although it is a 
well established measure, the mutual information does 
not distinguish between quantum and classical aspects 
of correlation, which is an important issue of quantum 
information theory. In recent years, the study of mea- 
sures that can differentiate among these two distinct as- 
pects of correlation has been intensified (T6l - [l8j , specially 
after the discovery that mixed unentaiigled states can 
also have non-classical correlation [l^, [TJ- The more 
widespread measure of total quantum correlation in a bi- 
partite quantum system is the quantum discord, which 
was proposed by OUivier and Zurek [l7| . A related quan- 
tity concerning classical correlations was also proposed 
by Henderson and Vedral [l^. A possible way to deter- 
mine the existence of correlated initial states was sug- 



gested by E.-M. Laine et al [20j and it is based on the 
trace distance, a measure of distinguishability between 
two quantum states [2l|. They have shown that if the 



trace distance increases over its initial value, the system- 
reservoir state was initially correlated. However, this ap- 
proach does not identify the nature of the initial system- 
environment correlations, that is, the quantum and clas- 
sical aspects. 

In this work we study the dynamics of a system con- 
sisting of a two-level atom (qubit) interacting with a zero- 
temperature environment without the Markov and Born 
approximations. We show that the system purity dy- 
namics can be used as a witness for initial correlations 
between system and environment. Furthermore, by intro- 
ducing a second atom (probe qubit) interacting with the 
same environment, we are able to determine the differ- 
ent initial system-environment correlations through the 
dynamics of the purity between the qubit and the probe 
qubit. We also study the dynamics of entropy and clas- 
sical and quantum correlations of the atomic system. In 
this case, we show that these dynamics can be useful to 
identify and to classify the initial correlations between 
system and environment. 



II. THEORETICAL MODEL 



In this section, we first describe theoretically some im- 
portant quantities that are employed to quantify quan- 
tum and classical correlations within quantum systems. 
These quantities are useful in our work because they 
exhibit signatures of initial system-environment corre- 
lations, as demonstrated by our results. Moreover, we 
present the model used to describe two atoms coupled 
to a common environment at zero temperature. The ap- 
proach to solve the dynamics of the system-environment 
is discussed as well. 
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A. Quantum and Classical Correlations 

Quantum Discord: Quantum correlations have played 
an important role in quantum information and communi- 
cation theory [15[ . If we consider for instance a bipartite 



quantum state described by the density matrix pab, the 
total correlations between the subsystems A and B can 
be computed by using the mutual information [l5j |: 



^ (pab) = S (pa) + S (pb) - S (pab) : 



(1) 



computational basis {|00) , |01) , |10) , |11)}. While all 
pure quantum correlated states are entangled states, 
there are quantum mixed states with null entanglement 
and a positive value of quantum discord. This result 
exemplifies that not all quantum correlations can be 
described by entanglement. 



B. Dissipative dynamics of two atoms coupled to a 
common environment 



where pi = Ttj {pij) and S{p) = — Tr (plog2 p) is the von- 
Neumann entropy. A measure of classical correlations 
present in the quantum state pab is defined as [Toj : 



J {pab) = S (pa) - S {pa\b) 



(2) 



where S {pa\b) = ^^^{Ah}J2bPbS{PA\b) is the quan- 
tum conditional entropy. The minimization is given over 
generalized measurements {Mf,} 15], with = Isi 

Mb > for all 6, and 



PAlb = 



{lA(g)Mb)pABilA(»Mb) 

Tr[{lA(E)Mb)pAB(lA<»Mb)] 



is the reduced density operator of A after obtaining the 
outcome b in B. For two-qubit systems, the minimization 
over generalized measurements is equivalent to a mini- 
mization over projective measurements [22i] . In this work, 
we use numerical minimization to compute the classical 
correlations. 

Using the results above, the quantum discord \17^ is 
defined as 



D {pab) = I {pab) ~ J {pab) ■ 



(3) 



and it assumes equal values, irrespective of whether the 
measurement is performed on the subsystem A or B, 
only when S (pa) = S {ps) [III 

Entanglement: A special class of quantum correlated 
states are the entangled states. A quantum bipartite 
state PAB is said entangled if and only if it cannot be 
written as a separable state pab = TlijPjPf ® ■> where 
J2jPj — 1 Hi- In this work we use the entanglement 
of formation (EoF) [23] as a measure of entanglement, 
which is defined as EoF{pab) = niin ^ -pjiSd^/ij)), 
where the minimization is over all ensembles of pure 
states {Pj,\4'j)} such that pab — J2jPj iV'i) ii^jl- For 
a pair of qubits, the EoF is a monotonically increasing 
function of the concurrence C [1^ 

EoF (pab) = -/(C)log2/(C)-(l-/(C))log2(l-/(C)), 

(4) 

where f{C) = (l + Vl-C^) /2 and C = 
max{0, Ai — A2 — A3 — A4} with Ai, A2, A3 and A4 
the square roots of the eigenvalues, in decreasing order, 
of the matrix R = p (ct^ (g) cr^) p* (cr^ (g) 0-^). Here p* 
denotes the complex conjugation of the matrix p in the 



In this subsection, we describe the dynamics of 
two non-interacting two-level atoms coupled to a zero- 
temperature common environment - the damped Jaynes- 
Cummings model. The environment is represented by 
a bath of harmonic oscillators with the Hamiltonian 
given by {h = 1) Hr = J2k^kalak, where al (ofc) is 
the creation (annihilation) operator for the field mode 
k with frequency uJk- Our system consists of two two- 
level atoms whose Hamiltonian (in computational basis) 
is Hq = Wo (c^+C- +cr^(T^), where ojq is the Bohr fre- 
quency and cr;^ and ai_ are the Pauli raising and lowering 
operators of the j-atom, respectively. The Hamiltonian 
of the system plus environment is = Hs + Hr + Hs r 
with 



HsR = (o"+ + cr+) ^ gkak + h.c, 



(5) 



where is the coupling constant. 

The dynamics can be solved using auxiliary variables 
defined from the properties of the S2ectral distribution 
through the pseudo-mode approach 26] . Here we assume 
that the two atoms interact resonantly with a Lorentzian 
structured reservoir, resulting in the following pseudo- 
mode master equation in the interaction picture |26| 



dp 
'dt 



- (2apa1" 



a'^ap 



pa'' a) 



(6) 



where p is the density operator for the two atoms and the 
pseudo-mode, a (a^) is annihilation (creation) operator 
of the pseudo-mode and V is the interaction Hamiltonian 
of the atomic system and the pseudo-mode. The spectral 
distribution associated with the pseudo-mode is 



J(w) 



TT (w-~wo)2 + r74' 



(7) 



where F is the pseudo-mode decay rate and fl is the 
coupling constant between the pseudo-mode and the two 
atoms. In this model, there are two regimes to be consid- 
ered [2^: (i) a. strong-coupling regime that occurs when 
F < 2r2, and (ii) a weak-coupling regime for F > 2U. 
As discussed in Ref. 26] , for strong-coupling regime, the 
atomic dynamics presents non-Markovian features and 
when the coupling becomes weak, the Markovian behav- 
ior is recovered. In all cases studied in this work, we 
assume that the total system has at most one excitation, 
so that the Hamiltonian V can be written as 1261 



V = ni 



h.c. 



(8) 
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This approach does not rely on either the Born or the 
Markov approximation and therefore it is possible to 
solve the dynamics without performing any further ap- 
proximation. 

III. RESULTS AND DISCUSSION 

One important issue in the theory of open quantum 
systems is the ability of determining the origin of the 
decoherence of a qubit H^l- For instance, if we start 
with the single qubit mixed state 

Ps{0) = a'\g){9\ + {l-a')\e}{e\, (9) 

where a G (0, 1), we are not able to determine how this 
state was initially prepared. For example, the qubit state 
described by Eq. ^ can be obtained from different situ- 
ations: 

a) the state was prepared without any correlation with 
the environment 

Pa = Ps(0)® |0)s(0|, (10) 

where \n) ^ means the state of the environment with n 
excitations; 

b) the system was classically correlated with the envi- 
ronment 

= a' 1.9) {g\mE + - «') Ie) {e\mE (0| ; (H) 

c) the system had quantum correlations and null en- 
tanglement with the environment 

Pc = a' \9) {9\mE i'f'l + i^ - Ie) (e|®|0)s (0| , (12) 

with |</))^ = (|0)^ + |1)^) 772; 

d) the system had non-null entanglement with the en- 
vironment 

P,(0) = |*(0))(vI'(0)|, (13) 

being |vI/(0)) =a|g)|l)^ + yT^|e)|0)^. 

In all states above, tracing over the environment vari- 
ables, we end up with the same density operator given in 
Eq. ^ . Then, as the information about the environment 
is not directly accessible, any measure carried out on the 
qubit state can not distinguish whether the system was 
initially correlated to the environment or not. Accord- 
ing to Ref. , the dynamics of the qubit can be used to 
identify the presence of initial correlations among system 
and environment. However, there is only the possibility 
of determining if the system was initially correlated to the 
environment or not. Based on this approach, we might 
say if the qubit initial state described by Eq. ^ was pre- 
pared in that way or it was a result from the interaction 
with the environment. 

In this work, we employ the dynamics of the system 
purity to witness the system-environment initial correla- 
tion. The purity P of a quantum state p is defined as 



P (p) = Tr (p^) and a quantum state is pure if and only 
if the density operator p is idempotent, that is, p^ = p. 

We begin our analysis by studying the dynamic behav- 
ior of the purity P of a qubit due to the interaction with 
the environment. We numerically calculate the dynamics 
of P for the initial states given in Eqs. p^ - ([T^ . for both 
strong and weak coupling regimes. In Fig. [U one can 
notice that the behavior of the dynamics of the system 
purity for the uncorrelated state is different from the cor- 
related ones. However, the dynamics of the entangled ini- 
tial state pd coincides with the classically correlated state 
Pb- These states differ from each other by two extra terms 
that appear in pa- For our type of interaction, these ex- 
tra terms evolve in time not affecting the dynamics of the 
reduced system, where the environment has been traced 
out. In this situation, the dynamics of the system pu- 
rity for both ph and pd becomes indistinguishable. Based 
on our results, we see that the dynamics of the system 
purity can be used as a measure to determine the system- 
environment initial correlation, however the kind of ini- 
tial system-reservoir correlation is not well determined. 
Moreover, we observe in Fig.[ljb) that the curves become 
indistinguishable when the system-environment coupling 
is weak, that is, the initial correlation does not interfere 
too much in the system dynamics, so the supposition of 
a non-correlated initial state is well justified in this case 
In addition, the purity asymptotically assumes the 
same value {P = 1), regardless the initial state. This 
is expected since we are assuming a zero-temperature 
reservoir and thus the asymptotic state of the qubit is 
the ground state \g). Although the Fig. [prefers to a par- 
ticular value of the parameter a — 0.5, our results were 
checked for different values of a between and 1 and the 
same characteristics were observed. 

By employing a second qubit (probe) interacting with 
the same environment and considering the initial states 
as p = \g) {g\ ® pi {i ~ a, b, c, d), we numerically evaluate 
the dynamics of the purity for the two qubits as shown in 
Fig. [21 For this situation, the dynamics of P have differ- 
ent behavior for each initial state. Therefore, we are able 
to determine the initial correlation between system and 
environment by observing the system purity dynamics. 
To understand this phenomenon, we need to analyze the 
role of a common environment. Although both qubits do 
not interact with each other directly, they are correlated 
because both qubits interact with the same reservoir. 
Such behavior has been observed in different physical 
systems. [1^ Moreover, the common reservoir produces 
correlations between the qubits and these correlations de- 
pend on the initial state of the system-reservoir. Such a 
dependency is more evident when the interaction between 
the qubits and the reservoir becomes stronger. Thus, the 
stronger the interaction, the better the distinction among 
the dynamics for each initial state. Besides, the initial 
correlation of the original qubit and the environment can 
be dynamically transferred to the second qubit in such 
a way that the purity of both qubits can produce a sig- 
nature of the kind of initial correlation among the first 
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FIG. 1: (color online) Time evolution of the one-qubit purity 
for (a) strong-coupling F = $1 and (b) weak-coupling F = 5il 
considering the initial states pa (solid line), pt = pd (dotted 
line) and pc (dashed line). We fixed = 0.5. 



FIG. 2: (color online) Time evolution of the two-qubit purity 
P for (a) strong-coupling F = and (b) weak-coupling F = 
50, considering the initial states pa (solid line), pt (dotted 
line), pc (dashed line) and pd (dash-dot line). We fixed = 
0.5. 



qubit and the environment. In other words, the two extra 
terms of the initial state pd induce correlations that are 
captured by the whole system (two atoms) purity. Again, 
in the asymptotic regime, the values of these quantities 
do not depend on the initial state. In Fig. [3J we plot 
the (a) mutual information, (b) classical correlation, (c) 
entanglement of formation, and (d) quantum discord for 
two-qubits as a function of dimensionless time fit for the 
same initial states of Figs. [T] and [H As already seen in the 
purity dynamics, each initial state leads to a particular 
dynamic behavior of the correlation measurements. An- 
other interesting aspect about these quantities is their 
usefulness to improve the identification of a particular 
state. For instance, the dynamics of the purity for the 
states Pc and pd are very similar (see Fig. [2]), which makes 
difficult the determination of the initial correlated state. 
On the other hand, the dynamics of the other quantities 
shown in Fig. [3] for the same states pc and pd are very 
distinct, thereby improving the ability of specifying the 
initial correlated state. 

To determine if there is or not initial correlation be- 
tween the qubit and its environment, it would be enough 
to compare the dynamics of the qubit purity obtained 
experimentally with the theoretical result. If it is nec- 
essary to obtain a more precise information about the 
type of the initial correlation, we must employ a probe 
qubit interacting with the same environment and then 
compare the dynamics of one of those physical quantities 
discussed above with the theoretical results. In order to 
reduce errors between experimental outcomes and the- 




OX Qi 

FIG. 3: (color online) Dynamics of the (a) mutual informa- 
tion, (b) classical correlation, (c) entanglement of formation, 
(d) quantum discord for two-qubits and for pa (solid line), 
Pb (dotted line), pc (dashed line) and pd (dash-dot line) with 
= 0.5 e F = Q. 



oretical results, the purity would be the best choice at 
first glance, since it is just obtained by matrix multiplica- 
tion (Tr [p^ (i)] ) instead of diagonalization or maximiza- 
tion processes, which must be performed to calculate the 
other physical quantities. If the results obtained from the 
purity are not decisive, one should perform an analysis 
of the other physical quantities to be able to determine 
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the initial correlated state. Thus, by employing such a 
procedure described above, one is able to infer if there 
is initial system-environment correlation and which is its 
nature. 

It is interesting to remark what would occur if other 
kinds of interactions in the Hamiltonian were considered. 
For example, an important class of such Hamiltonians is 
found when the system is coupled to a dephasing envi- 
ronment, where such a coupling is described by the Pauli 
matrix cr^. This kind of interaction is very common to 
model th e dy namics of solid states systems, e.g., quan- 
tum dots [29| and Josephson junctions [s^. If we consider 
this type of interaction together with the initial condi- 
tions treated in this work, we are not able to identify the 
initial correlation between the system and the environ- 
ment. This shortcoming occurs because the initial states 
of Eqs. (fT0|) - (IT3|) are written as an ensemble of eigenvec- 
tors of <7z, thus no difference on the dynamics will be 
observed for this type of interaction. However, if other 
initial conditions were considered, e.g., initial states rep- 
resented by an ensemble of eigenvectors of a^, the same 
methodology could have been applied. 

IV. CONCLUSION 

In summary, we studied a way to determine the exis- 
tence of initial correlation between a two-level atom and a 
Lorentzian structured reservoir at zero temperature. We 



found that the dynamics of quantum and classical corre- 
lations have signatures of the initial system-environment 
correlations. Moreover, we showed that it is possible to 
identify the distinct initial correlated states by consid- 
ering a probe qubit interacting with the same reservoir. 
We verified that, in our case, the purity might be the 
best choice to witness the initial system-environment cor- 
relations because this quantity reduces errors when the 
experimental outcomes are compared to the theoretical 
results, although the purity is not always decisive to de- 
termine the nature of the initial correlations. In such 
a situation, the dynamics of classical or quantum cor- 
relations measures for the bipartite system (two atoms) 
can be used to distinguish among different initial states. 
Finally, as we have only presented a numerical study, 
we hope that our work stimulates new researches in this 
area, specially in order to obtain analytical results that 
probes the efficiency of purity and the measures of cor- 
relations as witness to initial system-environment corre- 
lations for more general states and Hamiltonians. We 
also observed that both the purity and the measures of 
correlations become less efficient when the interaction be- 
tween the qubits and the reservoir is weakened, so that it 
is necessary to find a more precise witness to the initial 
qubit-reservoir correlation in the weak coupling regime. 

This work is supported by FAPESP and CNPq through 
the National Institute for Science and Technology of 
Quantum Information (INCT-IQ). 
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